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Abstract 

The Quantum Ergodic Conjecture equates the Wigner function for a typical eigenstate of a clas¬ 
sically chaotic Hamiltonian with a delta-function on the energy shell. This ensures the evaluation of 
classical ergodic expectations of simple observables, in agreement with Shnirelman’s theorem, but 
this putative Wigner function violates several important requirements. Consequently, we transfer 
the conjecture to the Fourier transform of the Wigner function, that is, the chord function. We 
show that all the relevant consequences of the usual conjecture require only information contained 
within a small (Planck) volume around the origin of the phase space of chords: translations in 
ordinary phase space. Loci of complete orthogonality between a given eigenstate and its nearby 
translation are quite elusive for the Wigner function, but our local conjecture stipulates that their 
pattern should be universal for ergodic eigenstates of the same Hamiltonian lying within a clas¬ 
sically narrow energy range. Our findings are supported by numerical evidence in a Hamiltonian 
exhibiting soft chaos. Heavily scarred eigenstates are remarkable counter-examples of the ergodic 
universal pattern. 
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I. INTRODUCTION 


The time average of any observable evolved by an ergodic classical Hamiltonian is deter¬ 
mined by the delta distribution over the energy shell. The Quantum Ergodic Conjecture 
(QEC) of Voros |I] and Berry [2] states that the same distribution can be reinterpreted as 
the semiclassical limit of the Wigner function P], for an eigenstate of the corresponding 
quantum Hamiltonian that has the same eigenenergy. 

As a consequence of this conjecture, the expectation value of a quantum observable is 
close to a classical average in phase space. Indeed, Shnirelman’s theorem HHE] conhrms 
asymptotically the conjecture for almost all eigenstates, albeit with severe restrictions for 
the allowed classes both of the observables and of the ergodic systems, (see e.g. Refs. HIE] 
and references therein). It should be noted that the QEC is stronger than Shnirelman’s 
theorem, with consequences that are much more subtle than the gross correspondence of 
expectations to their classical values. Most notable, so far, is the deduction by Berry |2| of 
local statistical correlations of wave functions, that are directly applied to the interpretation 
of experiments in quantum dots (see Ref. |9| and references therein). The ever increasing 
experimental refinement in manipulating quantum states may allow one to access other 
delicate properties, such as interference phenomena, for eigenstates of classically ergodic 
Hamiltonians. Thus, a more refined quantum ergodic conjecture with further predictions 
within a more precise framework may become a practical necessity. 

There are several features of a quantum state, specially a pure one, that are not correctly 
portrayed by the Wigner function proposed by the QEC. They were recently reviewed in 
[To], where it was shown that the operator that it represents is not even positive. Here we 
focus on the behavior of quantum states subject to a phase space translation. Its effect on 
the QEC Wigner function varies radically, depending on the choice of equivalent formulas 
that can be employed to describe it, pointing to yet another inadequacy of the QEC. 

Thus, we employ the QEC Wigner function as a stepping stone to formulate an alterna¬ 
tive Local Quantum Ergodic Conjecture (LQEC) for the Fourier transform of the Wigner 
function, in a neighborhood of the origin. LQEC encapsulates all desirable classical prop¬ 
erties of the standard QEC for the Wigner function. No assumption needs to be made 
outside this region, so there are no pitfalls. Furthermore, it predicts an universal pattern 
for the overlap of eigenstates, within a classically small energy window of a given ergodic 
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Hamiltonian, with their translation by a small displacement. 

This paper is organized as follows: In section after some preliminary dehnitions, the 


LQEC is presented. In section III, under a restriction on the Hamiltonian, we show that the 
LQEC implies the same results as the QEC for the wave function correlations for ergodic 


states. In section |IV[ the universal pattern of blind spots for ergodic states is presented, 
and tested in the Nelson Hamiltonian. Finally, we discuss our results in section |Vj The 
semiclassical theory near to the origin is summarized in the Appendix]^ Intermediary steps 
for the wave function correlations are contained in the Appendix 


II. LQEC AND SMALL CHORDS 

Denoting by x = (p, q) G a phase-space point, the Wigner function |3], IT(x), 
and its Fourier transform (known as the chord function [H]), x(^), are both complete 
representations of the density operator, p [m. They are conveniently dehned in terms of 
the unitary translation operator 




exp 




( 1 ) 


and in terms of the Fourier transform of T^, the reflection operator such that 

1 


X(^) = 


(flnh) 


D 


and 


IF(x) = 


flh) 


D 


tr P), 


tr {R^ p). 


( 2 ) 


( 3 ) 


According to the QEC, the Wigner function that represents a typical F-energy eigenstate 
of a classically chaotic Hamiltonian can be approximated by the classical Liouville probability 
density for the F-energy shell of the Hamiltonian with Weyl symbol hf(x) |12j . 

5{E - H{^)) 


WM = 


( 4 ) 


/ dx (5(F-id(x))’ 

where 6 denotes the Dirac delta-function. The assumption is that the classical Hamiltonian 
is ergodic, so that the time average of any classical observable, Ac{'x.), evolved by id(x) = 
H{p, q) according to Hamilton’s equations, coincides with the classical phase space average. 


(A) = / dx IT(x) A(x) , 


( 5 ) 
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for the classical probability density Q. The conjecture is then that one can reinterpret 
this same ergodic distribution, ITV(x), as the approximate quantum Wigner function that 
represents the £'-eigenstate of the corresponding quantum Hamiltonian. 

Let us focus on the behavior of quantum states subject to a translation. The overlap of 
a pure state with its translation is 

l(V'lV'^)r = l(^lAl^)l^ = (27rh)^y dx lT(x + ^) lT(x), (6) 

so that in the case of Q this would decrease monotonically to zero, until there is no more 
overlap between the energy shell and its translation. However, this scenario is contradicted 
by an alternative expression for the expectation of the translation of a pure state, that is, 

^ m = j w(x) = {i-m. (7) 

Indeed, insertion of Q into this formula reveals a hne structure of oscillations very close 
to the maximum at the origin of the phase space of chords, including zeroes, i.e. zero 
overlap between the state and its translation. 

This apparent paradox is resolved by shifting of the QEC conjecture for the Wigner 
function to the LQEC conjecture for the chord function: 

Xe(«) ^ j W,(x) (8) 

Here, both limits, ^ 0 and h —)■ 0, should be interpreted in a heuristic sense. This 

is familiar in the use of the semiclassical limit, in which good approximations may follow 
from the extrapolation to hnite values of Planck’s constant of results derived in the strict 
asymptotic limit, h —)■ 0. It will be shown that such extrapolations also hold for the short 
chord limit: The region in which the LQEC is presumed to hold is a Weyl volume of order 
surrounding the origin ^ = 0. 

The expression of the chord function in terms of its moments that is presented in the 
Appendix A implies that all moments calculated from the LQEC are classical, as are the 
expectation values of any observable that is polynomial in the components of position and 
momentum. Furthermore, the hrst few moments allow for a rough estimate of the location 
of the zeroes of x(^) close to the origin, generalizing the discussion in Ref. [13]. It is shown 
that there is a nodal surface for the imaginary part of x(^) that contains the origin and 
is locally parallel to (x). On the other hand, the lowest quadratic moments, forming the 
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covariance matrix [T3], K, lead to an approximate nodal surface for the real part of the 
chord function in the form of the little ellipsoid, namely 

((^Axf)=^ J^KJ^ = 2h^ (9) 

where J is the standard symplectic matrix and T stands for the transpose. These codimension- 
1 nodal surfaces intersect in codimension-2 surfaces, which have zero dimension in any two 
dimensional section of the chord function, so they generalize the blind spots, studied in Ref. 
[I3] for the special case where D = 1. 


III. WAVE FUNCTION CORRELATIONS 


So far, no restriction has been made on the momentum dependence of the Hamiltonian, 
but the second order Schrodinger equation follows from quadratic momenta, i.e. the Hamil¬ 
tonian separates into a quadratic momentum term and a potential energy, na 

ff(x)=p2 + V(q). (10) 


Then, the approximate chord function near the origin is 




(27rh) 


D 


dq 





6{E-p^-V{q)) 

J dp 6{E - p2 - U(q))’ 


( 11 ) 


so that, integrating over the surface of the p-sphere with radius r£;(q) 
obtains 


where we dehne 


Xe{$,) 


{27rh)^ He 


dq e 


fiSv 


^r?(|^q|;q,^), 


q, E) 


T{D/2) 


Jd_i [rEiqMci\/h] 
[rE{q)\Q/‘2h]^~^ 


^E - H(q), one 


( 12 ) 


(13) 


and Xe(^) = 0 otherwise. Here, J,y(x) denotes the Bessel function of order v [12]. Then, the 
U-volume for the region V (q) < E is just 


Hq = y dq En(0; q, E), (14) 

since En(0;q,E) = 1. For D = 2, one has simply Fod^ql; q, .E) = Jo(’^£;(q)|^q|/^)- 

The chord function, within the plane = 0, is merely the Fourier transform of the 
characteristic function for the region U(q) < E, divided by the H-volume of this region. 
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In the special case of a billiard, for which Id(q) = 0 within its bonndaries, there is no q- 


dependence in F, so that it factors out of (12). Then, the dependence of the chord function 
on resides only in F, divided by the H-volume of the billiard, and this is multiplied by the 
^p-dependent Fourier transform of the characteristic function for the interior of the billiard. 

The local chord function (12) is sufficient for the derivation of Berry’s (local) wave func¬ 
tion correlations [2], dehned as 




dq 


(q-Q) 


(27rA)^/2 ^ (q +V 2 |^)(^|q-V 2 ). 


(15) 


ArA(Q) 

where A/a(Q) is a normalization factor such that C'a(0, Q) = 1. This correlation function 
can be obtained as a Fourier transform of the chord function in a small neighborhood of the 
origin, as shown in the Appendix [Bj 

For the case of ergodic eigenstates, one inserts the approximate chord function of the 


LQEC ([^ into (B2), leading to 

Ca(«„Q) ^ '' 


A'a(Q) J (2irfi)D 


HE - g(x)) 


J dx 5iE - H{x)) 




(16) 


Then, assuming the Hamiltonian in the form (10), so that the chord function is given by 


(12), the correlation is given directly by (B2) as 


(Q-q) 


C'A(^q,Q) 


/dqe |;q,E) |;Q,E) 


(q-q)2 

f dq e FD{0]q,E) 


FniO] Q, E) 


(17) 


where the last approximation holds under the assumption that A is classically small, so that 
one can neglect the variation in H(q ~ Q). This is certainly the case for the local classical 
approximation of the chord function within a Planck volume. No assumption about the 
outlying regions of the chord function is required, so as to re-derive the isotropic wave 
function correlations of [2]. In the limit of a billiard, i.e H(q) = 0, the averaging region is 
not constrained, so that local correlations can be identihed with the global correlation for 
the wave function. 


IV. THE UNIVERSAL LOCAL CHORD FUNCTION 

Our general assumption is that the chord function for an energy eigenstate within a Planck 
volume of the origin is determined by corresponding classical structures at its eigenenergy. 
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In the case of an integrable Hamiltonian of a system with two or more degrees of freedom, 
neighboring states correspond to quantized tori that need not he close to each other within 
their neighboring energy shells. Hence, the local chord functions for neighboring states 
may be quite diverse. This will also hold for the pattern of their zeroes, that is, loci of 
possible translations of each eigenstate, where it becomes orthogonal to itself. The LQEC 
predicts an alternative scenario for the displacements of ergodic eigenstates: They all share 
a common classical structure, the full energy shell itself, and hence, for each Hamiltonian, 
the loci of zeroes of their local chord functions is universal, within a classically narrow 
energy range. Sweeping through two dimensional sections of the codimension-2 manifolds 
of zeroes, these zeroes arise generically as blind spots, just as for a single degree of freedom; 
the pattern of these blind spots is then predicted to be universal for neighboring ergodic 
states. Furthermore, the pattern of nodal lines for the real and the imaginary part of the 
local chord function should be reasonably approximated by Eq. Q. 

Any averaging over the chord function, or the Fourier transform of the Husimi function, 
masks the pattern of its zeroes. Therefore, the best way to visualize them is by examining 
two dimensional sections of the phase space of chords, that is, the locus of all possible 
translations. Then, the zeroes are typically isolated points of the section, blind spots just as 
in the case D = 1, lying at the intersection of the nodal lines of the real and the imaginary 
parts of the chord function [T3] . 

In the case that the potential energy, and hence the Hamiltonian, has a center of (reflec¬ 
tion) symmetry, the chord function is real, and hence the zeroes in any section are identihed 
with the nodal lines of the real chord function. Moreover, the chord function is then re¬ 
lated by a trivial re-scaling to the Wigner function itself irzi. Considering that a reflection 
symmetry is in no way incompatible with classical chaos, this is a dramatic falsihcation of 
a fully detailed account of the QEC for the Wigner function, which lacks the highest peak 
of any symmetric Wigner function at the reflection center (the origin), together with the 
surrounding oscillations. In contrast, this is precisely the region to which the LQEC applies 
and is here studied numerically. 

Note that Eq. (|^ can be evaluated numerically by Monte Carlo integration as 

1 ^ 

^ i=i 

where are random phase-space points sampled according to HQ(x). In practice, the 
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sampling is performed using a sufficiently narrow Gaussian about the energy-shell, namely 


WM ^ 




a 


exp 


(g(x) - Ef 

2cr2 


(a —)■ 0 ^ 


(19) 


We will consider the classical Nelson Hamiltonian is: 

-f^(Xl,X2) = - - -^ Y 


2 \ 2 
Ti 


( 20 ) 


and its quantum counterpart results from the replacement: Xi = {qi,pi) —t and 


X 2 = (^ 2 ,^ 2 ) —t (^ 2 ,^ 2 ) in (20). The classical dynamics of this system is mixed for all 
energies, but we consider eigenstates \'ipn) whose eigenenergies correspond to an essentially 
chaotic classical behavior (see Ref. [TH] and references therein). 

The exact quantum calculations are performed using the imaginary-time nonuniform 
mesh method (ITNUMM) proposed by Hernando et ah [I9]. This method is specially 
designed for the calculation high-excited states and it has a very easy implementation. 

For simplicity, we restrict to the plane of the momenta ^p, dehned by = 0. In a small 
(Planck) region about the origin. Fig. shows the LQEC and the exact chord function for 
the 109*^ excited state of the Nelson Hamiltonian in the plane of momenta. We observe 
that the LQEC furnishes an adequate description of the chord function in this small region 
about the origin. The blind spots arise as intersections of the nodal surfaces of the real and 
imaginary parts of x(^) fh® appropriate section. A comparison between of the exact and 
the LQEC nodal lines of the chord function for 109**^ state are shown in Figs. I^b) and (c). 
The closest blind spots to the origin can be roughly estimated by the intersection of the 
the short-chord ellipsoid ([^ and the plane ^ A (x) = 0. In the ^p-plane, the latter is the 
horizontal axis, whereas as the former is depicted by the green line in Fig. |^b). Figures]^ 
show the exact quantum chord functions for two consecutive states, the 109*'^ and the 110**^, 
respectively. Here a universal structure in a neighborhood of the origin is observed; however, 
it is not present in regions far from the origin. The chord function in ([^ decay exponentially 
in the outer region. This pattern depends only on the shape of the energy-shell. 

For scarred states, the LQEC and the exact quantum differs substantially. Fig. |^a) 
shows the 39*^ excited state of the Nelson potential and the respective scarring classical 
orbit. We observe that the nodal pattern of this scarred state differs noticeably from that of 
typical ergodic states at nearly the same energy, so that the nearest blind spots are displaced, 
as show in Figs. |^b) and (c). 








Real 


Imaginary 



FIG. 1. Exact quantum (upper) and semiclassical LQEC (lower) chord functions for the 109 


ith 


state of the Nelson Hamiltonian (20) in the ^p-plane. Here h= 1. 


Exact (0.6895) LQEC (0.6310) - Short-chord (0.4792) 



EIG. 2. 109*^ excited eigenstate of the Nelson Hamiltonian (20). (a) Exact wave function, (b) 

Real and (c) imaginary nodal lines of the chord function in the ^p-plane, evaluated by several 
methods: exact quantum (red dashed), semiclassical LQEC (gray solid) and short-chord ellipse 
(green solid). In parenthesis is the estimation of the first blind spot along the horizontal axis for 
each method. Here h= 1. 


9 


































































Real 


Imaginary 



FIG. 3. Exact chord functions for the 109**^ and 110*^ excited state of the Nelson potential in the 
^p-plane. Here h= 1. 


Exact (1.085) LQEC (0.7016) -Short-chord (0.8278) 



FIG. 4. Scarred state of the Nelson potential for the 39*^ state: (a) The density probability in 
position space {qi,q 2 )- Red line shows the classical periodic orbit, (b) Real and (c) imaginary 
nodal lines of the chord function: exact quantum (red dashed), semiclassical LQEC (gray solid) 
and short-chord ellipse (green solid) 
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V. DISCUSSION 


The subtle nature of the transition from the QEC to the LQEC may beneht from an 
alternative view point: The hard core of quantum ergodicity is Shnirelman’s theorem and so 
concerns expectations of mechanical observables, whose Weyl symbols are smooth functions. 
For this end the QEC is perfectly adequate, even though one cannot consider Eq. Q to 
be a proper Wigner function. By including the real and the imaginary part of the chord 


function for small chords as the expectations of the Hermitian operators (see Eq. (A5) in the 
Appendix]^ as legitimate observables for the QEC, we substantially increase the predictive 
value of our ergodic conjecture, the LQEC, by encompassing quantum features that are 
more rehned than simple expectation values. This small region of the chord function allows 
for the re-derivation of wave function correlations that were previously known and predicts 
blind spots for small chords that are quite problematic for the QEC. 

Perhaps, the greatest advantage is that the LQEC is not burdened by a full prediction 
that cannot be dismembered within the Wigner picture. On the contrary, it is just within 
a tiny region of the chord phase space that the ergodic chord function is determined, so 
that one avoids the violation of general rules. By the identihcation of the LQEC with the 
small chord approximation [TB], in the case of systems with a single degree of freedom, one 
obtains a clear notion of the complete irrelevance of the extrapolation of the LQEC beyond 
the Planck volume that it legitimately describes. 

The value of the outer region of the chord function is to reveal the purely quantum 
features that distinguish quantum states. In the case of neighboring ergodic eigenstates, the 
expression 

|(V’n|V'n')P = (2vrh)^ [ Xn{^) Xn'i^T = 0, (21) 


reveals that the similarity of both chord functions close to the origin must be compensated 
by their differences further off. These differences, which can be discerned in Fig. then 
highlight the gap in our knowledge of ergodic eigenstates. 
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Appendix A: Moments of position and momentnm 


The dehnition of the chord fnnction Q allows ns to calcnlate the distribntion moments, 
or statistical moments, of p and q in the form of the derivatives of the chord fnnction, i.e. 
explicitly 


(p") = Tr p^p = {-ihy 


QUx 




^=0 


(Al) 


and 


(q’^) = Tr = {iKy 


QUx 




(A2) 


^=0 


Conversely, if we know all the moments, then we know the chord fnnction, becanse the 
expansion in a Taylor series of the chord fnnction is 




Mb 

(ih)^ 



{M (q’-^p*)) 


n—k 
P ’ 


where, for conjngate components {qj,Pj) 






(A3) 


(A4) 


and Pnk stands for all possible permntations of prodncts of g” and p^. (No permntation is 
needed for commnting components.) 

Thus, knowledge of the chord function in an inhnitesimal neighborhood of the origin is 
equivalent to information on the expectation for all observables in the form of polynomials 
in position and momenta. In principle, the whole chord function would be determined by 
all statistical moments combined with analiticity, this is no available in practice, due to the 
very oscillatory nature of highly excited eigenstates. 

Recall that the chord function ([^ is in general complex. Because it represents a Hermitian 
operator, then x(—^) = x(^)*; where the asterisk denotes complex conjugation. The real 
and imaginary parts of the chord function are related to the cosine and the sine transforms 
of the Wigner function, respectively. By dehning the Hermitian operators = (T^ + T_^)/2, 
and sj = (T^ — T_^)/2i, one can see that 


^cx(^) oc {yjlc^iy) and 3mx(^) oc (AS) 


12 








Appendix B: Small chords and wave function correlations 


Inserting the inverse chord transform m, 


(q +W2|^)(^|q-V2) = 






into the local wave function correlation (15), one obtains 

1 


Ca(«<„Q) = 


Vi(Q) 

Note that the normalization constant can be expressed as 

Vi(Q) = /<«p x(0.«p) 


(Bl) 


(B2) 


(B3) 


Let us then choose the smoothing parameter as A i.e. the width of a coherent 

state, so that the Gaussian integration window in Eq. (15) has a volume Within 


this classically small volume, the wave function for a highly excited eigenstate necessarily 
involves an average over many oscillations of wave length h. Then, the same Gaussian width 


is obtained for the smoothing of the chord function in (B2), so that the local dehnition of 


the wave function correlation only requires the information of the chord function within a 
Planck volume surrounding the origin. It should be noted that the chord function of a state 
that is translated in phase space by a chord r/. 


\',(«) = x(«) 


(B4) 


leads to the interpretation of (B2) as the Gaussian smoothed projection of the chord function 


for the eigenstate translated by r/ = (Q,0). 

Gonsider the effect of a further smooth projection, but now onto the P-axis with the 
specihc choice A = 


^ A/'a(Q) GA(^q,Q) = 


[ttH) 


D 


(nh) 


^-$y2h-i. 


D 


=w(x). 


(B5) 


Thus, we identify the Husimi function [201 HI], "H, evaluated at the phase space point, 
X = (Q,P), that may be dehned as the Gaussian coarse-graining of the Wigner function. 
Surprisingly, the product of a Gaussian amplitude on the chord function does not diminish 
its information content, since the Husimi function is also a complete representation of a 
quantum state. Here we hnd that exact knowledge of the local correlation for all positions. 
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Q also contains, in principle, complete information about the state. In the case of the 
Husimi function, information is delicately encoded in its analytic properties and may be 
hard to extract. 
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